HEP-01-13 



A Class of Exactly- Solvable Eigenvalue Problems 



Carl M. BenderH and Qinghai Wangi 

Department of Physics, Washington University, St. Louis, MO 63130, USA 

(February 7, 2008) 



Abstract 

The class of differential-equation eigenvalue problems —y"{x) + x'^^^'^y{x) = 
Ey{x) {N = — 1, 0, 1, 2, 3, . . .) on the interval — oo < x < co can be solved 
in closed form for all the eigenvalues E and the corresponding eigenfunctions 
y{x). The eigenvalues are all integers and the eigenfunctions are all confluent 
hypergeometric functions. The eigenfunctions can be rewritten as products of 
polynomials and functions that decay exponentially as x — > ±00. For odd 
the polynomials that are obtained in this way are new and interesting classes 
of orthogonal polynomials. For example, when = 1, the eigenfunctions are 
orthogonal polynomials in multiplying Airy functions of x"^. The properties 
of the polynomials for all N are described in detail. 
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I. INTRODUCTION 



In this paper we consider the class of differential-equation eigenvalue problems 

- y"{x) + = x"" Ey{x) (iV = -1, 0, 1, 2,3,.. .) (1.1) 

on the interval — oo < x < oo. The eigenfunction y{x) is required to obey the boundary 
conditions that y{x) vanish exponentially rapidly as x ^ ±cxd. For each integer > —1, it 
is possible to solve these eigenvalue problems in closed form. The eigenvalues are all integers 
and the associated eigenfunctions are all confluent hypergeometric functions. Furthermore, 
all the eigenfunctions for each value of N can be written as the product of a polynomial 
and a given function that vanishes exponentially for large The classes of polynomials 
that are obtained in this way are orthogonal and for odd are apparently new and have 
interesting mathematical properties. 



The eigenvalue problem ( |1 . 1| ) discussed in this paper arises in many contexts. In classical 
physics a perturbative technique called boundary- layer theory has been developed to find 
approximate solutions to boundary-value problems of the form 

ew"{x) + a{x)w\x) + h{x)w{x) w{-l) = A, w{l) = B, (1.2) 

where e is treated as a small parameter. Problems of this sort appear in the study of fiuid- 
fiow problems in various geometries. Perturbative treatments of this equation are usually 
quite straightforward |T|. However, there is a particularly difficult special case of ( |1.2| ) that 
may occur when there is a resonant internal boundary layer. Suppose that a(0) = 0, so that 
there is an internal boundary layer at x = 0. Suppose further that near x = 0, a{x) ~ ax 
and h{x) ~ (3. Then, near x = the differential equation in ( |1.2|) is approximated by 

ew"{x) + axw\x) + I3w{x) = 0. (1.3) 

The Gaussian change of variables w{x) = e~°'^^^^^'^''y{x) converts ( |1.3| ) to Schrodinger form: 



Apart from a scaling, this equation is the A^ = case of ( |1.1D . It describes the quantum 
harmonic oscillator, and eigenvalues occur when the parameters a and (3 satisfy 

(3={n + l)a (n = 0,1,2,3,...). (1.5) 

When the parameters a and (3 are related in this fashion, the internal boundary layer is 
said to be resonant. Unlike conventional boundary layers, a resonant boundary layer is not 
narrow; that is, its thickness is not small when e is small. As a result it is particularly 
difficult to treat the resonant case using ordinary boundary-layer methods |l|]. 

More generally, if a{x) ~ ax^^^ and h{x) ~ I3x^ when x is near 0, we obtain the 
differential equation in ( p.2|) . This equation is then converted to the differential equation 
in ( |1 . 1| ) by the exponential change of variables in ( ^.1[ ). Thus, the eigenvalue problem in 
( |1.1| ) characterizes the general resonant case in the theory of internal boundary layers. The 
special case of the harmonic oscillator discussed above occurs when A^ = 0. 
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The eigenvalue problem also arises in the context of supersymmetric quantum mechanics 
and quasi-exactly solvable models. It appears, for example, in the recent work of Voros 
and Dorey et a/ 0. The same differential equation was also considered by Znojil 0] 
but with boundary conditions imposed on a semi-infinite interval. The quantum problem in 
( p.. ID may be thought of as a peculiar inverse approach to quasi-exact solvability. Ordinarily, 
in this field one tries to construct potentials for which a finite number of eigenvalues of the 
spectrum can be calculated exactly and in closed form, while the remaining part of the 
spectrum remains analytically intractable. The problem in (|1 . 1| ) is to construct potentials 
V{x) for which there is an eigenvalue that is exactly zero. The zero eigenvalue may or 
may not be the ground-state energy of the potential V{x) = x^^"^^ — x^E that has been 
constructed. In this paper we will see that the case of odd-integer is much more interesting 
than the even-A^ case. 

The eigenvalue problem in (|1.1|) is especially interesting because, as we show in Sec. ||, 
leading-order WKB theory (physical optics) gives the exact spectrum E for all odd N and 
almost the exact answer for even N. 

This paper is organized as follows. In Sec. |T| we give the exact solution to the eigen- 
value problem in (|1 . 1|) . We show that the eigenvalues E are all integers for each value of 
N = —1, 0, 1, 2, 3, . . . and that the corresponding eigenfunctions are all confluent hyperge- 
ometric functions. In Sec. ^ we examine the eigenfunctions for even-integer N. For this 
case the eigenspectrum is positive and the nth eigenfunction has definite parity. The nth 
eigenfunction has the form of a polynomial of degree n and of argument x^"*"^ multiplied by 
the exponential exp(— x^"^^), which decays as x ±oo. The polynomials are generalized 
Laguerre polynomials. The polynomials for even n form an orthogonal set and the polyno- 
mials for odd n form a different orthogonal set. In Sec. IV, we study the eigenfunctions for 
odd-integer N. For this case the spectrum of eigenvalues E ranges from — oo to oo and the 
nth eigenfunction does not exhibit definite parity. For each the nth eigenfunction has 
the general form xA^(x^+i)P„(x^+2) + A^(x^+i)Q„(x^+2). Here, P„(z) and Q„(z) are 
polynomials of degree n that satisfy the same recursion relation but have different initial 
conditions. The functions Ai\i{z) are independent of n and are generalized Airy functions 
that obey the differential equation A'l^{z) = z'^^^'^'^^^ Aj\i(z). When A^ = 1, the function Ai{z) 
is the conventional Airy function Ai{z). In Sec. we consider the special cases A^ = — 1 and 
A^ = 1 (the Airy case). 

We emphasize that the eigenvalues that are obtained in this paper are not the energies 
of a conventional Schrodinger equation. However, the eigenfunctions that are obtained 
might well be useful for solving some conventional Schrodinger equations. For example, it 
might be useful to express the eigenfunctions of the pure anharmonic oscillator problem, 
—y"{x) + x^?/(x) = Ey{x) as linear combinations of the eigenfunctions found in this paper. 
The work done in this paper suggests that it may well be advantageous to expand the 
function y{x) as P(x)Ai(x^) + Q(x)Ai'(x^), where P and Q are series in powers of x. 



II. EXACT SOLUTION OF THE EIGENVALUE PROBLEM 

We solve for the eigenvalues of the differential equation ( |1.1| ) by converting it to a con- 
fluent hypergeometric equation and then imposing the boundary conditions. We begin by 
making the substitution 
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1 „iV+2 

y{x) = e ^+'2 w{x). (2.1) 
The function w{x) then satisfies the differential equation 

w"{x) - 2x^+^w'{x) + (3x^w{x) = 0, (2.2) 

where 

f] = E~N-l. (2.3) 

There are now two cases to consider, even and odd. Suppose first that N is even. 
There are two hnearly independent solutions to the differential equation ( |2.2| ): 

^(a;) = ^F, (-—JL^, 1 —■ .2.4) 

^ ^ ^ ^ 2(A^ + 2)' A^ + 2' A^ + 2 / ^ ^ 

and 

When the first parameter of the confiuent hypergeometric function is a negative integer, 
its Taylor series truncates to a Laguerre polynomial: 

,F,{-n,c;t)=n\-p^Lt'\t). (2.6) 
i c — n 



Thus, for the solution in ( |2.4| ) we obtain 

/9n 



2{N + 2) 



n (n = 0,1,2,3,...), (2.7) 



or 



E„ = 2n(Ar + 2) + Ar + 1 (n = 0,1,2,3,...). (2.8) 
The corresponding eigenfunctions have even parity. For the solution in (|2.5|) we obtain 

- -J— + , , = n (n = 0,1,2,3,...), (2.9) 

A^ + 2 2(iV + 2) ^ y j 

or 

Er, = 2n{N + 2) + N + 3 (n = 0, 1, 2, 3, . . .). (2.10) 

The corresponding eigenfunctions have odd parity. 

When is odd, the solution to the differential equation (^I^) that is well behaved as 
X ^ oo is a particular linear combination of the solutions in ( p.4| ) and ( p.5| ) known as a 
confluent hypergeometric function of the second kind: 
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[2.11] 



Note that as x — > +00, the function y{x) in (2J.) vanishes exponentially. However, as 
X — > — 00, the exponential factor in ( |2.1| ) grows and we have 



X 



^N+2 



N + 2 



1 



2{JV+2) 



N+2 



r 1- 



1 



Af+2 



X 



1 + 



Af+2 2(Ar+2) 
1 \' 



r 1- 



1 



+ 



/3 



N+2 ' 2{Ar+2) 



r 



2{N+2) 



r(i + 



/3 



2(Af+2) 



\X 



N+2 



(2.12) 



Because exp (— ]v+2-^^^^) S^^i'^s exponentially as x —00, the only way to satisfy the 
boundary condition as x — >■ — cxo is for the expression in square brackets to vanish. The 
expression in square brackets simplifies to 



sm 



\N^ 



+2 2{N+2) 



IT 



sm 



2{N+2) 



sm 



{n+2) 



(2.13) 



Hence, 



1 



N+2 N+2 
and thus we obtain the eigenvalues 

E„ = (2n + l)(Ar + 2) 



+ 2n (n = 0,±l,±2,±3,...). 



{n = 0,±1,±2,±3,...). 



(2.14) 



(2.15) 



It is interesting that a leading-order WKB analysis (the physical optics approximation) 
of (|1.1|) gives the exact eigenvalues for odd and almost the exact eigenvalues when N 
is even. Consider the following two-turning-point time-independent Schrodinger equation 
boundary-value problem 



y"ix) + Qix)y{x) = 0, i/(±oo) = 0. 



(2.16) 



Ordinarily, Q{x) = V{x) — E, where ^(x) is the potential and E is the energy. In the 
physical-optics approximation, the condition for a solution to this problem to exist is 



j dx \J —Q{x) 



n + -]tc 



(n = 0,l,2,3,...). 



(2.17) 



where the turning points A and B satisfy Q{A) = Q{B) = 0. If we apply the quantization 
condition (|2.17|) to (|1 . 1|) , where Q{x) = x^^^^ — Ex^ , we obtain 
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dx ^eJ^'^^i^^ = + -K (n = 0,1,2,3,...), (2.18) 

where we assume without loss of generahty that E is positive. The turning point B satisfies 
]^N+2 _ This integral can be evaluated exactly as a beta function and we obtain 

E = (2n+ l)(A^ + 2). (2.19) 

This is precisely the result in ( ^.151 ) for odd. Also, it is nearly the results for even N in 
( p.8|) and ( |2.1CI|) , which can be combined to read 

E = (2n+ l)(A^ + 2) ± 1. (2.20) 

A striking property of the eigenvalue problem ( p..l[ ) is that for even A^ the eigenvalues are 
positive but for odd A^ the eigenvalues are both positive and negative. This is reminiscent 
of the difference between the Klein-Gordon equation for bosons, which has positive-energy 
states only, and the Dirac equation for fermions, which has positive-energy states (electrons) 
and negative-energy states (positrons or holes). 

In the next two sections we describe the two cases A^ odd and A^ even in greater depth. 



We consider the simpler case of even A^ in Sec. |Tl| and turn to the more interesting case of 
odd A^ in Sees. fV] and 0. 



III. EIGENVALUE PROBLEM FOR EVEN N 

When A^ is even, all the eigenvalues are positive and the eigenfunctions have either even 
or odd parity. The even-parity eigenfunctions have the form 

y2n{x) = e-i^^""^Ll-^) (ivT2'^''"'') = °' ^' ^' ^' • • ^^-^^ 

and the corresponding eigenvalues are 

E2n = 2n{N + 2) + AT + 1, (3.2) 
where L^") is the generalized Laguerre polynomial. The odd-parity eigenfunctions are 

y^n+iix) = e-i^-""'x lI^^ (n = 0, 1, 2,3,.. .) (3.3) 

and the corresponding eigenvalues are 

E2n+i = 2n{N + 2) + N + 3. (3.4) 

Note that the eigenfunctions have the form of a decaying exponential multiplying a poly- 
nomial. For the even-parity solutions we can write the polynomial as a monicQ polynomial 
Pn{z) in the variable z = 4x^"'"^: 



^The term monic means that the coefficient of the highest power in the polynomial is 1. 
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Pr^iz) = {-irn\ [2{N + 2)]« U ^-^> (^^^^^J (3.5) 
and for the odd-parity solutions we have 

qr^iz) = (-l)'^n! [2{N + 2)]" li^^ (2(7VT2) ) " ^^"^^ 
These polynomials satisfy very similar recurrence relations 

Pn+i{z) = [z- 2{N + 2)(2n + 1) + 2\pn{z) - 2{N + 2)n[2(iV + 2)n - 2]pn-i{z), 

qn+i{z) ^[z- 2{N + 2){2n + 1) - 2]q^{z) - 2{N + 2)n[2{N + 2)n + 2]g„_i(^), (3.7) 

where the initial conditions are 

PoW = l, px{z)^z-2N -2, 

go(^) = l, gi(^) =^-27V-6. (3.8) 
The polynomials Pn{z) and qn{z) also obey similar differential equations 
2{N + 2)zpl{z) + (27V + 2 - zy^{z) + np„(z) = 0, 

2{N + 2)z?:(^) + (27V + 6 - ^)?;(z) + nq^{z) = 0, (3.9) 

and differential relations 

2(iV + 2)zplSz) + Pn+i(^) + [2(iV + 2)(n + 1) - 2 - z\pn{z) = 0, 

2(Ar + 2)zq'^{z) + g„+i(z) + [2{N + 2)(n + 1) + 2 - z]qn{z) = 0. (3.10) 

The generating functions for these polynomials are also quite similar: 

CO / -i\nj.n ^ J 

Gp{z,t) = E , Pni^) = [l-2(iV + 2)t]-^+We™^, 

n=0 

G,{z, t) = Y^ ^ \ qn{z) = [1 - 2(iV + 2)t]-'-—^ eTO^. (3.11) 

n=0 ^■ 

The polynomials Pniz) and 5n(-2^) a-re separately orthogonal: 

1" dzwp{z)pm{z)pn{z) = [2(7V + 2)]2"n! ^^^^^^^ 5„.,n, 

° ^ JV+2 j 

J dzw,{z)qUz)qr,{z) = [2{N + 2)]'^n\^ 5m,n, (3.12) 

where the weight functions Wp{z) and Wq{z) are given by 
Wp(z) = e"5(i^2;~^[2(A^ + 2)]^"^ 



Wg(z) = e 5(j^^ivT2[2(Ar + 2)]~^"^— ^ (3.13) 



The moments of these weight functions are 

I V 



/ \ r°° T (n + l + 



a^f) ^ I dzw,{z)z- = [2{N + 2)]"^^^-^^. (3.14) 



The (divergent) power series constructed from these moments have particularly simple con- 
tinued fraction expansions in which the continued-fraction coefficients are all integers: 

00 -1 

V a^P^t" = = 

„t^o " 1 [2(iV + 2) - 2]t 



[2(iV + 2)]t 



[4(A^ + 2) - 2]t 



_ [6(A^ + 2) -2]t 

m+2)i 



^ = [2(iV + 2) + 2]t ■ ^^-^^^ 



n=0 



[2(Ar + 2)]t 



[A{N + 2) + 2]t 



[4(Ar + 2)]t 



[6(iV + 2) + 2]t 
^ [6(iV + 2)]t 



We illustrate these general results for the two special cases = and N = 2. 
Special Case = 0: The Harmonic Oscillator. For this case the eigenvalues in ( p.2|) 
are E = 1, 3, 5, 7, ... and the polynomials Pn{z) and xqn{z) in (|3.5|) and (|3.6| ) coalesce to 
become the standard Hermite polynomials H„(x): 

Pn{z) = R2n{x) = (- l)'^22"n! L^-i/^) (x^) , 
xq^iz) = H2„+i(x) = (-l)"22"+in!xL(i/2)(x^). (3.16) 



Special Case N = 2. For this case the eigenvalues in ( ^.2[ ) are = 3, 5, 11, 13, 19, 21, 
and the first few monic polynomials p„(x) and g„(x) in (|3.5| ) and (|3l6| ) are 

Poiz) = 1, 

Pi(z) = z-6, 

P2{z) = 2^-28^ + 84, 

P3(z) = - 66z^ + 924z - 1848, 

p^(z) = z^- I20z^ + 3960^2 - 36960Z + 55440, 

qo{z) = 1, 
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qi{z) = z -10, 

q^lz) = - 36z + 180, 

q^(z) = z^ - 78z^ + 1404z - 4680, 

q^(z) = z^ - 136z^ + 5304^2 - 63648^ + 159120. (3.17) 



IV. EIGENVALUE PROBLEM FOR ODD N 

For this case the eigenvalues are 

En = {2n + l){N + 2) (n = 0,±l,±2,±3,...) (4.1) 

and the corresponding eigenfunctions yn{x), which do not have definite parity, are confiuent 
hypergeometric functions of the second kind: 

yjx) = U (-n - — i -, 1 -^x^'-^A . (4.2) 

^ V 2(A^ + 2)' A^ + 2'A^ + 2 / ^ ' 

Note that the boundary condition as a; +oo is aheady satisfied and the quantization 
comes from requiring that y{x) — > as a; — > —oo. 

For each the eigenfunctions yn{x) in ( [4. 2] ) can be expressed in terms of what we will 
call generalized Airy functions An{x) combined with polynomials Pn and Qn as follows:^ 

JV+l 

y-n-i{x) = 2 ^^xAn 



and 

yn{x) = y-n-i{-x) (n > 0). (4.4) 

We define the generalized Airy function A]y{x) as the solution to the differential equation 

A'l^{x) = x^An{x) (4.5) 

that decays exponentially as a; ^ +oo. Note that when = 1, the function Ai{x) is just 
the conventional Airy function Ai(x). We can express An{x) in terms of the associated 
Bessel function K,,{z) as follows: 

An(x) = —UiN + l)]2WT2)J^^K iv+i f^^^x^V (4.6) 



^To avoid confusion, for odd N we use upper-case notation P„ and Q„ to represent the polynomials; 
we use lower-case notation p„ and qn to represent the polynomials associated with even A'^. 



iV+l 



2 2{iV+2) X 



N+1 



P„(4x 



N+2\ 



N 



N+l 



2 2(JV+2) X 



N+l 
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With this choice the function An{x) is normahzed so that 

f°° 1/A^+1\ / A^+1\ N+7 1 3 

dxA^ix) ^ -r l^-^j r (^ ^^^^ 2-TO(7V+ l)w(7V + 2)-OT. (4.7) 

Note that this reduces to the standard result /g°° dx Ai[x) = | when = 1. 

The polynomials Pn{z) and where z = 4x^~^'^, both satisfy the same recursion 

relation 

Pn+-,iz) = [z + 2{N + 2){2n + 1)]P„(^) - [2{N + 2)n - 1][2{N + 2)n + l]Pn-i{z), 
Qn+i{z) ^[z + 2{N + 2){2n + l)]Qn{z) - [2{N + 2)n - 1\{2{N + 2)n + l]Qn-i(^), (4.8) 

but have different initial conditions 

Po(-2) = l, Pi(^) = ^ + 2Ar + 5, 

QoW = l, Qi(^) =^ + 2Ar + 3. (4.9) 

The polynomials Pniz) and Quiz) satisfy coupled second-order differential equations 

4(iV + 2)zP::{z) + 4{N + ?,)P'^{z) + 2zQ'^{z) + Qr,{z) = (2n + l)P„(z), 

4(7V + 2)zQ"^{z) + 4(iV + \)Q'^{z) + 2;sP;(;2) + P^{z) = (2n + l)gn(;2), (4.10) 

and coupled differential relations 

4(iV + 2);^P^(^) = 2P„+i(z) - zQr^iz) - [A{N + 2)(n + 1) + 2 + z]P„(z), 

4(7V + 2)^Q;(^) = 2Qn+i{^) - zPn{z) - [4(7V + 2)(n + 1) - 2 + z]Qn{z). (4.11) 

The generating functions for these polynomials are 



oo j.n 

Gp{z,t)^Y.-M^) 



n=0 ^■ 



\Fl[l+ O.J, ox, 1 



[1 - 2(7V + 2)i]-^-TO ' ' ^ 2(iV+2)^ ^ ^ iV+2. 2(iV+2)[l-2(iV+2).]y ^ ^^^^3) 

iPl (^1 + 2(Ar+2) ' ■'■ Ar+2' 2(iV+2), 



and 



00 jn 

GQ(^,i) = ^-g„(^) 

n=0 ^• 



= fl - 2(7V + 2)il-^+TO 5^+^)' iv+2>2(iv+2)[i-2(iv+2)^]y_ .^_^3. 
'■ ^ ^ ' F M 1 1 _ 1 . z \ ^ ' 

1-^1 V 2(Ar+2) ' Ar+2' 2(Ar+2) j 

The polynomials Pn(-z) and Qj^z^ obey identical-looking orthogonality and normalization 
conditions 
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/oo 1 TT / 1 \ 

/oo 
-oo 



sin ,rTT^^[2{N + 2)]2"+ir (n + 1 



TT 2(7V + 2) 
xr ( n+ 1 + 



2(7V + 2)^ 



2(7V + 2)J 



(4.14) 



The weight functions Wp{x) and VFq(x) are real and positive and are expressible as principal- 
part integrals: 



Wp{x) 

Wgix) 



dsV 



dt 

dsV 

-OO J —oo t — S 



dsV / 



In 



In 



-t 



2(A^ + 2)7r 



e 4^^^ K iv+3 



-t 



2(Ar + 2)7r 

or, in terms of the generaUzed Airy functions A]^{x), 

fX foo I 1 t 

Wp{x)= dsV In^ -V27r(-t)5(i^e ^^^A 

J — OO J — OO t — 5 I 

/X roo (j^f I 1 t 

dsV In ^ V27r(-t)^(^^e~2(iv+2)^^ 
-oo J —oo t — S I 



5p?+5y \^4(A^ + 2)^ 



(4.15) 



N+l 

-t \ 



.4(A^ + 1) 



N+l ■ 
N+2 



(4.16) 



where V indicates principal-part integration and the integral is performed on the sheet for 
which — 1 = e^*'^. 

The moments of the weight functions Wp{x) and Wq{x) are given by 



/■GO roo 

4)=/ dxx^'Wpix) and a^^^ = dxx^'Wqix). 
J— 00 J— 00 



(4.17) 



The divergent power series constructed from these moments have remarkably simple 
continued-fraction expansions in which the continued-fraction coefficients are all integers: 

00 1 



n=0 I 



[2(iV + 2) + l]t 



[2(iV + 2)- l]t 



1 - 



[4(Ar + 2) -f l]t 



1 - 



[4(Ar + 2) - l]t 

[6(A^ + 2) + l]^ 
^ [6{N + 2)-l]t 
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n=0 



[2(A^ + 2)-l]t 



(4.18) 



[2(Ar + 2) + l]t 



[4(iV + 2) - l]t 



[4(iV + 2) + l]t 
[6(A^ + 2) - l]t 
[Q{N + 2) + l]t 



V. TWO SPECIAL CASES OF THE ODD-iV EIGENVALUE PROBLEM 

In this section we consider two interesting special cases of the odd-A^ eigenvalue problem; 
namely, = — 1 and = 1. 

Special case N = —1. For this case Eq. ([4.5|) is of course not valid. However, the formula 
for the eigenvalues E in ( |4.1| ) is still valid and E = ±1, ±3, ±5, ±7, . . .. The eigenfunctions 
yn{x) in (|4.2|) are now Bateman functions kE„{x): 

2 /•'^/2 

Unix) = ksAx) = — dd cos(x tan^^ — EnO). (5.1) 
7r JO 

Apart from an overall multiplicative constant, the solution can be written in terms of asso- 
ciated Bessel functions combined with polynomials. For negative eigenvalues 

y-n-i{x) = xK^{x)Pn{Ax) + xK'^{x)Qn{Ax) (n > 0), (5.2) 

and for positive eigenvalues 

yn{x) = y-n-i{-x) (n > 0). (5.3) 

Note that the eigenfunctions are finite at the origin but that there is a branch cut. For 
definiteness, we take the branch cut to run up the positive imaginary-x axis. 
In terms of the variable z = Ax the first few polynomials Pn{z) and Qn{z) are 

^o(^) = 1, 

Pi(z) = ^ + 3, 

P2(^) = ;2' + 9^ + 15, 

P^[z) = z^ + 19^2 + 90^ + 105, 

P^[z) = z^ + 33;z^ + 321^2 + 1050^ + 945, 

P^[z) = z^ + hlz^ + 852^3 + 5631^2 + 14175z + 10395, 

Qo{z) = 1, 

Qi{z) = z + l, 

Q2{z) = z^ + 7z + ?,, 

Q^(z) = z^ + I7z'^ + 58Z + 15, 

Q^(z) = z^ + 31/ + 261/ + 582z + 105, 

Q^[z) = / + 49/ + 756/ + 4209/ + 6927^ + 945. (5.4) 
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The polynomials Pn{x) and Qn{x) satisfy the recursion relations in ( [4.8| ) with = — 1, 
the coupled second-order differential equations in ( 4.10| ) with = —1, and the coupled 
differential relations in ( [4.11| ) with N = —1. The generating functions for Pn{z) and Qn{z) 
are expressed in terms of Bateman functions k^: 

oo j.n p4(l-2t) h „ ( 5 1 

Gp(.,t) ^ E -Mz) = (1 - 2t)-^'^ - 



n=0 



t 



GQ(^,t)^^-Q„(^) = (l-2t) 



A;_3(2/4) 

p4{l-2t) , ( 2 1 

1/2^ '^-H4(l-2t)J 



-1/ 



n=0 



n 



fc_i(2/4) 



(5.5) 



The polynomials Pn{z) and Qn{z) obey identical- looking orthogonality and normalization 
conditions 

rfa;iyp(a;)P„(a;)P„(x) = / dxWQ{x)Qm{x)Qnix) = (2n - 1)!! (2n + 1)!! (5.6) 

-oo J— oo 

where (—1)!! = 1. Note that the weight functions Wp{x) and Wq{x) are real and positive 
and are expressible as principal-part integrals: 



Wp{x) = r dsV r -^In .rie-'/'K,{~t/4) 
J-oo J-oo t — s V 27r 



Wq{x) 



dsV 



°o (it 



t- s 



In 




e-*/^iro(-t/4) 



(5.7) 



The moments of the weight functions Wp{x) and Wq{x) give rise to the following lovely 
continued-fraction expansions: 



n=0 



3t 



n=0 



5t 



3t 



1 - 



7t 



5t 



1 - 



3t 



3t 



5t 



5t 



It 



(5.8) 



Special case = 1 (Airy functions). For this case the eigenvalues E in ( [4.1D are 

= ±3, ±9, ± 15, ± 21, . . . and the eigenfunctions yn{x) in ( [4.2|) are written in terms 
of Airy functions combined with polynomials. For negative eigenvalues we have 
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y-n-i{x) = 2-^/3 ^ Ai(2-2/3x2)P„(4x=') +Ai'(2-2/3x2)g„(4x=^) {n > 0), (5.9) 

and for positive eigenvalues we have 

yn{x) = y.n-i{-x) (n>0). (5.10) 

The polynomials Pn and Qn are functions of the variable z = Ax^. The first few such 
polynomials are 

^'o(^) = 1, 

P,iz)=z + 7, 

P2{z) =^2 + 25^ + 91, 

p^[z) =z^ + 55^2 + 698^ + 1729, 

P^[z) =z^ + 97^^ + 2685^^ + 22970z + 43225, 

Qo{z) = 1, 

Qi{z) = z + b, 

Q2{z) =^2 + 232 + 55, 

Q^(z) =z^ + 53z^ + 602z + 935, 

Q^(z) =z^ + 95z^ + 2505^2 + 18790^ + 21505. (5.11) 

The polynomials Pn{x) and Q„(x) satisfy the recursion relations (|4.8| ), the coupled 
second-order differential equations ( |4.10| ), and the coupled differential relations ( ^.11| ) with 



= 1. The generating functions Gp{z,t) and GQ{z,t) and the integral representations of 
the weight functions Wp and Wq are obtained by setting = 1 in ( [4. 121 ), ( [4. 131 ), ( [4.15[ ), 



and (^4.16|) . The polynomials Pn{z) and Qn{z) obey identical-looking orthogonality and 
normalization conditions 

J^^dzWp{z)PUz)Pn{z)= J^^dzWQ{z)QUz)Qn{z) = h'^r (n + T (u + ^-^ 6^,n 

(5.12) 

and the moments of the weight functions Wp{x) and Wq{x) have the following continued- 
fraction expansions: 

oo -■ 



7t 

n=0 I 



1 ^ 

^ . 13t 



lit 
~T9t 
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oo 



E = • (5-13) 



5t 

n=0 



7t 



lit 



^ 13t 



17t 



19t 



There is an interesting connection between the moments and the combinatorial 
numbers 6*2^, which represent the sum of the symmetry numbers of the 2n- vertex connected 
vacuum graphs in a 0^ quantum field theory p: 

qW) = 6nCg4". (5.14) 
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